ANALYTIC CONTINUATION OF EIGENVALUES OF THE LAME 

OPERATOR 



KOUICHI TAKEMURA 

Abstract. Eigenvalues of the Lame operator are studied as complex- analytic func- 
tions in period r of an elliptic function. We investigate the branching of eigenvalues 
numerically and clarify the relationship between the branching of eigenvalues and 
the convergent radius of a perturbation series. 



1. Introduction 

The Lame equation is an ordinary differential equation given by 

d 2 



dx 2 



+ n(n + l)p(x))f(x) = Ef(x), (1.1) 



where p(x) is the Weierstrass p-function which is doubly-periodic with a basic period 
of (1,t), n G Z>i and E is a constant. In [8, §23] and §15] this equation is 
discussed in detail. 

To analyze the spectral of Eq. (jl.l|) . we can choose boundary conditions in various 
ways. One is to impose a non-trivial periodic or anti-periodic solution to Eq. (|l.ip . 
Then, the set of eigenvalues E is discrete and the periodic or the anti-periodic solution 
is called the Lame function or the singly-periodic Lame function. Another is to impose 
a non-trivial doubly-periodic solution to Eq. (jl.l|) . In this case the set of eigenvalues 
E is finite and the doubly-periodic solution is called the Lame polynomial. When we 
change the variable by z = p(x), the doubly periodic function is essentially expressed 
as a polynomial in z. Related to quantum mechanics we can choose a boundary 
condition to have a non-trivial square-integrable solution on the interval (0, 1) to 
Eq. (|l.l|) . We remark that the eigenvalue E with each boundary condition depends 
on the period r. 

In this paper we investigate how the eigenvalues of Lame functions depend on r. 
In particular we consider branching of the eigenvalues as a complex-analytic function 
in r for the case n — 1. 

Set q = exp(7rv / — It). It is shown in 6J that eigenvalues never stick together if 
q G R and — 1 < q < 1. Therefore if q G R and — 1 < q < 1 then there is no branching 
of the eigenvalue E function in q (or r). 

Also note that we can calculate eigenvalues of Lame functions as power series in q 
by considering perturbation from the trigonometric model (the case q = 0) as written 
in jHj. It is proved in jB] that the convergence radius is not zero. If the convergence 
radius is 1, the eigenvalue is analytic in r on the upper half plane, but it is observed 
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numerically that the convergence radii of some eigenvalues are not 1 (see section |3J). 
Hence there must exist a singularity on the convergence circle. 

On the other hand it is known that for the Lame equation with n G Z>x or more 
generally for the Heun equation with integer coupling constants, the global mon- 
odromy is expressed by a hyperelliptic integral [Zj. As an application we obtain a 
condition for q that causes a branching of eigenvalues of the Lame function (see |Zj or 
section 0] in this paper). By thorough calculation, we obtain numerically some values 
q which produce branching. 

Finally, we find that the absolute value of the branching point calculated by in- 
vestigating the hyperelliptic integral nearly coincides with the convergence radius 
calculated by perturbation expansion. In other words we obtain a compatibility be- 
tween the global monodromy written as a hyperelliptic integral and the perturbation 
expansion through the branching point. 

This paper is organized as follows. In section El we review several choices for setting 
the boundary conditions for the Lame operator and observe their relationship. In 
section El we explain results on perturbation and the convergence radius. In section 
0]we consider the global monodromy and search for branching points numerically. In 
section |3] we discuss the compatibility between perturbation and branching points. 
In the appendix, several propositions are proved and definitions and properties of 
elliptic functions are provided. 

Throughout this paper, we assume that n is a positive integer, and we use the 
conventions that f(x) is periodic f(x + 1) = f(x), f(x) is anti-periodic f(x + 
1) = —f(x) and f(x) is doubly-periodic <^> ((/(x + 1) = f(x) or —f(x)) and (f(x + 
r) = f( x ) or -f(x))). 

2. Boundary value problems of the Lame operator 
We consider boundary value problems of the Lame operator H, where 

H = -^+n(n + l)p(x). (2.1) 

Let <Ti n t(H) be the set of eigenvalues of H whose eigenvector is square- integrable 
on the interval (0, 1), i.e. 

a mt (H) = {E\3f(x) G L 2 ((0, 1)) \ {0}, Hf(x) = Ef(x)}. (2.2) 
Let <Jd(H) be the set of eigenvalues of H whose eigenvector is doubly-periodic, i.e. 
a d (H) = (2.3) 
{E\3f(x) + s.t. Hf(x) = Ef{x), f{x + 1) = ±f(x), f{x + r) = ±f(x)}, 

Note that the doubly-periodic eigenvector is simply the Lame polynomial. It is known 
[HI that #a d {H) = 2n + 1. 
Let a s (H) be the set of eigenvalues of H whose eigenvector is singly-periodic. Set 

a p (H) = {E\3f{x) + s.t. Hf{x) = Ef(x), f(x + 1) = f(x)}, (2.4) 
a ap (H) = {E\3f(x) jt s.t. Hf{x) = Ef(x), f(x + 1) = -f(x)}. (2.5) 

Then a s (H) = cr p (H) \[a ap {H). On the sets ai nt (H), ad{H) and <J S (H) we have 
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Proposition 2.1. (i) For r G R + y 7 — 1R>0; we have 

a int {H)Ua d {H) = a s {H). (2.6) 
(ii) Assume that q = exp^y 7 — It) G R and < |g| < 1. Then 

o- mt (H)U°d(H)=o- s (H), (2.7) 
i.e., a int (H) U cr rf (if) = a s (H) and a int (H) n cr d (i?) = </>. 

We prove this proposition in the appendix. Note that, if q is not real, then 
the proposition o- int (H) fl o~d(H) = might be false. In fact, if n = 1 and q = 
V 7 — 1(-3281 . . . ), then it seems that —e\ G ai nt (H) fl o~d(H) (see Proposition 14.21 and 
Table 3). 

Next, we briefly explain the relationship to the finite-gap potential. Let 

d 2 

I = -^ + n{n + l)p{x + r/2) (2.8) 

and <Tb(I) be the set such that 

E G o~b(I) -v^ Every solution to (I — E)f(x) = is bounded on x G R. 
Ince [2] established that, if q = exp^y 7 — It) G R, then 



\ a b (H) = (-oo, £ ) U (£ l5 £ 2 ) U • ■ ■ U (£ 2n _ 1; £ 2n ) (2.9) 



where <Jb(H) is the closure of the set cr&(i?) in C, £!$ G o~d(H) and E < Ei < ■ ■ ■ < 
E 2n . Hence there is a finite band structure on eigenvalues of unbounded eigenvectors. 
This is referred to as finite-band potential or finite-gap potential. 

3. Perturbation and convergence radius 

In this section we calculate eigenvalues of Lame functions as power series in q(= 
exp(7rv / — It)). For this purpose we consider perturbation from the trigonometric 
model. First we consider a trigonometric limit q — > (<=>■ r — > y/— loo) and later 
apply a method of perturbation from the trigonometric model. 

For the case q = the spectral problem becomes much simpler. Set 

d 2 , , 7T 2 , 
^ = --5-2+71(71 + 1)— . (3.1) 

dx z sin 7rx 

Then if — > ify — ^-n(n + 1) as g = exp(7rv— It) — ► 0. The operator ify is the 
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Hamiltonian of the Poschl- Teller system or the +i trigonometric Calogero-Moser- 
Sutherland system. Set 

$(x) = (shi7nr) n+1 , v m = c m C^ +1 (cos7rx)$(x), (m G Z> ), (3.2) 



where the function C^(z) = ^rj^x 2Fi(—m,m + 2v; v + ~; is the Gegenbauer 
polynomial of degree m and c m = ^J 2 + ^^^2n+2)^ n+1 ' > ■ Then 

H T v m = n 2 (m + n + l) 2 v m , (3.3) 
and (v m ,v m >) = 5 m , m ', where the inner product is defined by 



</, 5)= / f(x)g(x)dx. (3.4) 
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Set 

H = ^ / : R — > C measurable 



£ \f( x )\ 2dx < +°°' 

f(x) = f(x + 2) a.e. x, }, (3.5) 

f(x) = (-l) n+1 /(-x) a.e. x 



H+ = {/eH|/(i) = /(x + l) a.e. x}, 
H_ = {/eH|/(x) = -/(i + l) a.e. x}. 

Inner products on the Hilbert space H and its subspaces H + , H_ are given by (•,•)• 
Then we have H + J_ H_ and H = H + ©H_. The Hamiltonian H (see Eq. (|2.1Jl ) acts 
on a certain dense subspace of H (resp. H + , H_) and the space spanned by functions 
{v m \m G Z> } (resp. {v m \m G 2Z> }, {v m \m G 2Z> + 1}) is dense in H (resp. H + . 
H )■ 

Now we apply a method of perturbation and have an algorithm for obtaining eigen- 
values and eigenf unctions as formal power series of q. For details see 

Set q = exp(7r-\/— It). For the Lame operator (see Eq. (|2.1|) ). we adopt the notation 
H(q) instead of H . The operator H(q) admits the following expansion: 

2 00 

H{q) (= H) — H T — ^-n(n + 1) + ^ V 2k (x)q 2k , (3.6) 

k=l 

where Ht is the Hamiltonian of the trigonometric model and V2k(x) are functions in 
x which are determined by using Eq. (jB.7|) . 
Set 

7T 2 

E m = vr 2 (m + n + l) 2 n(n + l). (3.7) 

3 

Then v m is an eigenfunction of the operator H(0) with the eigenvalue E m . 

Based on the eigenvalues E m (m G Z> ) and the eigenfunctions v m of the operator 
H(0), we determine eigenvalues E m (q) = E m + Y^=i Em k ^q 2k and normalized eigen- 
functions v m (q) = v m + Y^h=i l>2m'£Z >0 c m% v m'ff k of th e operator H(q) as formal 
power series in q. In other words, we will find E m (q) and v m (q) that satisfy equations 

H{q)v m {q) = (H(0) + ££° =1 V 2k {x)q 2k )v m {q) = E m {q)v m {q), (3.8) 
(v m (q),v m (q)) = 1, 

as formal power series of q. 



First we calculate coefficients ^2 m , €Z &^Lv m > — V2k{x)v m (k G Z >0 , m G Z 



>o 



Next we compute Em^ and c^^, for > 1 and m,m' G Z> . By comparing coef- 
ficients of v m >q 2k in Eq. ()3.8jh we obtain recursive relations for Em and c£^,. For 

details see j^j. Note that, ilm — w! is odd, then = c^^, = 0. Convergence of 

the formal power series of eigenvalues in the variable q obtained by the algorithm of 
perturbation is shown in [0]. 

Proposition 3.1. [SJ Corollary 3.7] Let E m (q) (m G Z> ) (resp. v m (q)) be the formal 
eigenvalue (resp. eigenfunction) of the Hamiltonian H(q) defined by Eq. \3.tf\) . If \q\ 
is sufficiently small then the power series E m (q) converges and as an element in the 
Hilbert space H the power series v m (q) converges. 
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We show an expansion of the first few terms of the eigenvalue E m (q) and the ra- 
dius of convergence for the case n = 1 in Table 1. We calculate the expansion of 
E m (q) = E m + ^2 k Em k ^q 2k for more than 100 terms and approximate the abso- 
lute values of coefficients Em k ^ by ab 2k for some constants a and b which are deter- 
mined by the method of least squares. Then, the radius of convergence is inferred 
by liminffc-too 1/ {\Em k ^\/ a) 1 ^ 2k . The inferred radius of convergence and expansions 
of the first few terms of the eigenvalue E m (q) are calculated as follows: 



„2/lU | 80 n 2 I 1360 n 4 j 20800 n 6 j 195920 n 8 i 3174880 n 10 j 684960 n 12 j \ 1 7zL q 

n r 3 " 27 ^ 243 " 2187 " 19683 " 59049 " ' ' ' / 

2 / 4fct j 27g 2 I 198928 4 , 55403584 6 4307155408 " 8 i 2879355070048 10 j \ j7K~ 

I .-I. l&H M&JL 759.375, ^ ~T 341 71 875,. ,-f,_"t, .38443359375 ^ ' / " ' ^ y 



E (q) 



EM 
EM 



_2 I WO I 09^2 1 22Yy24S 4 1 i. , I > I 6 i 1654Y62851USS „8 _i_ \ 07c 

11 \ .3 Z .35 1 Z 42875 52521875 ^ Z 12867859375 ~l~ • • • J ( •) 

77; — o — : — 77Z — a — : — rrn — n — : — vrnx. — a — : — Trm — Tn — : — ttott^ — i 



7T 2 (f + 20q 2 + 65g 4 + + W + W u + ^q 12 + ^g 14 + . . . ) TS3^ 

tt 2 I + fg 2 + + ^-7" • • • ^m^'r + • • • ) -838 



Ei(q) 



EM 
EM 



7l 2 



1 27 " 486 " .34992 " 629856 " 7558272 ± ' ' ' I 

rggl I 82 2 I 50339 4 I 13640101 6 I 3872868499 "'8 3267409458867 'lO , ^ nnc 

L-3_ZL_5_z_ZLaotKLz_ZL 200000 g .32000000 ^ Z .32000000000 1 Z : : u I 



Table 1. Expansion of the first few terms and the inferred radius of convergence. 

We introduce propositions on the spectral of the Hamiltonian H on the Hilbert 
spaces for the case q 2 G E and |g| < 1. Let an(H) (resp. <7h + (H), cth-(H)) be the 
spectral of the operator H on the space H (resp. H + , H_). 

Proposition 3.2. (c.f. [6, Propositions 3.2, 3.5] J Let q 2 G IR and |g| < 1. TTie opera- 
tor H is essentially selfadjoint on the Hilbert space H (resp. H +; H_J. The spectrum 
o"a{H) (resp. o~n + {H), o~-h._(H)) contains only point spectra and it is discrete. 

Proposition 3.3. (c.f. [6, Theorem 3.6] j Let q 2 G R and |g| < 1. v4// eigenvalues of 
H on the space H can be represented as E m (q) (m G Z> ), which is real-holomorphic 
in q 2 G (—1, 1) and E m (0) = E m . The eigenfunction v m (q) of the eigenvalue E m (q) 
is holomorphic in q 2 G (—1, 1) as an element in L 2 -space, and the eigenvectors v m (q) 
(m G Z> ) form a complete orthonormal family on H. 

It is shown that, if g 2 G M, |g| < 1 and m G 2Z> (resp. m G 2Z> + 1), then the 
corresponding eigenvector v m (q) belongs to the space H + (resp. H_) and we have 

o- H (H) = {E m (q)\m e Z> } (3.9) 

= {E m (q)\m G 2Z> } 
a H _(#) = {£ m (g)|mG2Z> + l} 

Among the spaces <th(H), <th + (H), o~n_(H), a int (H), o~ p (H) and a ap (H), the fol- 
lowing relations are satisfied: 

Proposition 3.4. VFe have crn(H) = <T% n t{H), o~n + (H) = a int {H) n o~ p (H) and 
g-u.( h ) = a int (H) n a ap (H) 

Proof. It follows from the definition of H that, if f(x) G H, then the function f(x) 
is square-integrable on (0,1), i.e. crn(H) C cr int (H). Now we show o int (H) C 
an(H). Let G amt(H). Then there exists a non-zero function /(x) such that 
Hf(x) = Ef(x) and f Q \f(x)\ 2 dx < 00. The exponent of the differential equation 
(H-E)f(x) = at x = is {— n, n + 1}. Since the function /(a;) is square-integrable 
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and the equation (H — E)f(x) = is invariant under the transformation x <-» —x, 
the function f(x) is expanded as 

f(x) = x n+1 (c + c l x 2 + c 2 x i + ...) (c ^0) (3.10) 

and satisfies f(x) = (—l) n+1 f(—x). From the periodicity, the function f(x + l) is also 
an eigenfunction. The function f(x + 1) is written as a linear combination of f{x) 
and another linearly independent solution, and we have f(x + 1) = C f(x) for some 
0) because f(x) is locally square-integrable near x — 1. It follows immediately 
that f(-x - 1) = C~ x f{-x). We have C/(x) = /(x + 1) = (-l) n+1 /(-x - 1) = 
(-l) n+1 C- l f(-x) = C- l f(x). Hence C G {±1} and f(x + 2) = /(x). Therefore we 
have f(x)eH,Ee an(H) and o int (H) C (?h(H). 

Relations <th + = cr int (H) f]a p (H) and cr H (H) = a int (H) r\a ap (H) are obtained 
by considering periodicity □ 

It is shown that eigenvalues never stick together as in j^j. 

Proposition 3.5. (c.f. [HI Theorem 3.9]) Let E m (q) (m e Z> ) be the eigenvalues 
of H(q) defined in Proposition \3.^A If q 2 G R and \q\ < 1, then E m (q) ^ E m /(q) 
(m 7^ m'). In oi/ier words, eigenvalues never stick together under the condition 
q 2 G R and \q\ < 1. 

Proof. Assume that the proposition is wrong. Then there exists m and g such that 
E m {q) = E m+ i(q). Let /(x) and f(x) be the corresponding eigenfunctions. Then 
one of f(x) or f(x) is periodic and the other is anti-periodic. Hence f(x) and f(x) 
are linearly independent. Since there is no first differential term in H, we have 

(£f( x )) />) - = °- Hence (s/W) f( x ) - f( x )£f( x ) is a constant 

and it is non-zero by linear independence. It contradicts the periodicity of f(x) and 
f(x) and we obtain the proposition. □ 

Corollary 3.6. (c.f. Corollary 3.10],) If q 2 G R, \q\ < 1 and m < m! , then 
E m (q) < E m ,(q). 



4. MONODROMY AND BRANCHING POINTS 
We consider the monodromy of solutions of 

Hf(x) = Ef(x), H = -^ + 2p(x) (4.1) 

dx z 

for each E. Note that this is the case n — 1 in Eq. fjl.l}) . 

For the case n = 1, we have o~d(H) — {— e i> — e 2, — e 3 } and the corresponding 
doubly-periodic eigenfunctions are pi(x), p2(x), ps(x) (see Eq. (|B.5|) ). From the peri- 
odicity of pi(x) (i = 1, 2, 3) we have Cd(H) PI cr p (H) = {— ei} and crd(H) n a ap (H) = 
{-e 2 ,-e 3 }. 

We now consider the expression of solutions to Eq. (j4.1j) for each i£. The functions 
H(x, £7) and -P(-E") defined around Proposition lA.ll for the case n = 1 are calculated as 
E(x,E) = p(x)+E and P(E) = (E + ei )(E + e 2 )(E + e 3 ). Then the function A(x, E) 
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defined in Eq. flA.4|) is a solution to the differential equation (jl.l|) (see Proposition 
IA.2|) . and it is also expressed as 

A(x, E) = A^!e-<W, £ = -pfo), (4.2) 

for suitably chosen A (see jU §39] or [HI §23.7]), where cr(a;) is the Weierstrass sigma- 
function and ((x) is the Weierstrass zeta-function (see Appendix). Note that we 
can show directly that the function A(x,E) written as Eq. (J4.2j) satisfies Eq. (J4.1j) . It 
follows from Eq. (|4.2j) and Eq. ()B.3|) that the monodromy is described as 

A(x + 1, E) = A(x, E) exp (2^ - C(*o)) , (4.3) 

where 771 = C(V 2 )- Hence, if 2^it - C(*o) e vrv^lZ (resp. 2^ - C(*o) e 27r v /T TZ, 
2771*0 - C(*o) e 27r v /r TZ + vrv^T), then £ G a s (H) (resp. £? G 0- p (#), E G a ap {H)). 
It follows from Proposition IA. 31 that . if 2rjit — ((t ) G" it\/—YL, then £ G" ai nt (H). By 
Proposition 12. II and Proposition 13.41 if —1 < q(= exp(7Tv / — It)) < 1, then we have 

Oh(#) = ^(#) \ {- e i> - e 2, -e 3 }, (4-4) 
<r H+ (#) = \ {-ei}, *h-(H) = (Tap(H) \ {-e 2 , -e 3 }. 

The eigenvalue in cr p (H) is analytically continued in g (or r) as to preserve the 
property 

E = -p(t ), 2^*0 - C(t ) e 2^^^. (4.5) 
and the eigenvalue in a ap (H) is analytically continued in q (or r) as to preserve the 
property 

E = -p(to), 2^*0-^0) e27rV = TZ + 7rV=T. (4.6) 
It follows from the relation i£ = — p(to) and Eq. ()B.4|) that Eq. (|4.3|l is rewritten as 

A(x+l,E)=A(x,E)expl~f E -2 m d ~\ ^ 

V J - £l V /-(^ + e 1 )(E + e 2 )(E + e 3 ) / 

Hence we reproduce the monodromy formula in terms of (hyper)elliptic integral which 
was obtained in |7j. For analyticity of elements in c p (H) or a ap (H), we have 

Proposition 4.1. (c.f. Theorem 4.6 (ii)]) If the eigenvalue E satisfies Eg. (|^-5| ) 
or Eq. \4-oV , E — 2r\\ 7^ and E 7^ — e±, —e<i, —e^ at q = q*, then the eigenvalue E 
satisfying Eq. (|^.5| ) or Eq. \4-$ ^ s analytic in q around q = q*. 

Note that Proposition 14.11 is proved by applying the implicit function theorem as 
is done in Theorem 4.6 (ii)]. The following proposition describes the condition for 
q (or r) that the set o~d(H) H ai nt (H) is non-empty. 

Proposition 4.2. Under the assumption E G o~d(H) (i.e., E G {— e±, — 62, —£3}), the 
condition E G ai n t(H) is equivalent to the condition E — 2i]i = 0. 

Proof. It follows from the assumption that E = —Ci for some % G {1,2,3}. A solu- 
tion to Eq. (j4.1j) for E = — is written as pi(x), and another solution is written as 
pi(x) j(l/pi(xf )dx. By Eqs. (jB~3llB~6|) we have 

dx f dx f (p(x + ujj) — e,i)dx ((x + Ui) + eiX 

(4.8) 



Pi(x) 2 J p(x) - ej J (e* - e v ){ei - e,») (e* - e^e* - e v >) 



8 



KOUICHI TAKEMURA 



where G {1,2,3} with i' < i", i ^ i', and i ^ i". Set si(x) = pi(x) and 
S2(x) = pi(x)(((x + Ui) + tiX — T]i). Then they are a basis of solutions to Eq. (j4.1|) 
for E = —ei, and si(x) (resp. s 2 (x)) is odd (resp. even). Since Si(x) has a pole at 
x = and s 2 (x) is holomorphic at x = 0, square-integrable eigenfunction on (0, 1) is 
written as As2(a;) for some constant A. Since £2(2 + 1) cannot have a pole at x = 
for square-integrability and it is written as 

s 2 (x + 1) = + 1) (C(x + uJi + 1) + (x + l)e; - jfc) (4.9) 
= ± (s 2 (x) + (ei + 2r} 1 )p i (x)) 

for some sign ±, we have E — 2rji = (i.e., — — 2rji = 0). 

Conversely, if E — 2r\\ = and E = — then it follows from Eq. (j4.9|) that s 2 (a;) is 
perioic with a period 1 and it is holomorphic on M. Hence s 2 (x) is square-integrable 
on (0, 1), and we have E G ai nt (H). □ 

By Propositions 12. ll 14.11 and 14.21 it follows that if the eigenvalue E in u p {H) or 
a ap (H) has a branching at q, then we have E — 2rji = 0. Hence a necessary condition 
that the eigenvalue E in <J P {H) or a ap (H) has a branching is that g and to satisfy the 
following conditions: 

2 Vl = -p(t )(=E), (4.10) 

27/ito - C(*o) e vrv/^TZ. (4.11) 

We try to solve Eqs. ()4.10l 14. 1 If) numerically. First we fix the value q. We expand 771, 
p(t ) and C(*o) m <? according to Eq. (|B.7|) with approximately 100 terms, and solve 
Eq. (|4.10j) numerically by Newton's method and obtain t . We evaluate Eq. (|4.11[) 
using to and check whether it is satisfied or not. Note that the imaginary part of the 
value to should be taken to be small in order to exhibit good convergence. 

By investigating more than 1000 complex numbers which satisfy \q\ < .90, > 
and Qq > where Kg (resp. 5s q) is the real part (resp. the imaginary part) of the 
number q, we obtain numerically that the numbers in Table 2 may have branches 
(i.e. they satisfy Eq. ()4.10|) and Eq. (|4.11Jl ). Note that it seems some numbers do not 
generate branching. 



periodic 


q = .328106/, .258666 + .6974481, .510303 + .546057/ 

.746852 + .452463/, .224582 + .842777/, .552288 + .677536/ 
.314813 + .821858/, .686317 + .559106/ 


anti-periodic 


q = .281417 + .534362/, .655163 + .503275/, .264829 + .792687/ 
.535905 + .640487/, .807197 + .405705/ 



Table 2. Numbers which may have branches. 



Next we consider how to continue the eigenvalues analytically in q along a path. 
Let C be a path in the complex plane. The eigenvalue E is continued analytically in 
q along the path C by keeping the conditions 

E = -p(t ), (4.12) 

3m G Z, 2?7ito - C(*o) = m-n^f^l. (4.13) 

Note that the eigenvalue satisfying Eq. (j4.12j) and Eq. (j4.13|) for m G 2Z (resp. m G 
2Z + 1) is continued from the eigenvalue in H + (resp. H_). 
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We solve Eqs. (|4.121 14.13J1 for points which are selected appropriately on the path C 
and are connected by choosing close solutions. Note that for each E and q satisfying 
Eqs. (l4.12( 14.13(1 . solutions (to,m) may not be unique. Sometimes we need to change 
to another solution (t' ,m f ) to avoid the divergence of continued solutions in q. 

We continue the eigenvalue E analytically around the possible branches in Table 2. 
We obtain that the following numbers would not cause branching and they all would 
satisify 2r)i = — for some i G {1, 2, 3}: 



Q = 


.328106/ 


2r/i = 


-ei 


q = 


.281417 + 


.534362/ 


277! = 


-e 2 


Q = 


.510303 + .546057/ 


2rn = 


-ex 


Q = 


.655163 + 


.503275/ 


277! = 


-e 2 


Q = 


.746852 + .452463/ 


2th = 


-ex 


Q = 


.264829 + 


.792687/ 


277i = 


-e 3 










Q = 


.807197 + 


.405705/ 


2??i = 


-e 2 



Table 3. Numbers that do not cause branching. 

For these cases, it is inferred from Proposition 14. 21 that one of the eigenvalues E m (q) 
(m G Z> ) meets with an eigevalue with doubly-periodic eigenfunction (i.e. — ex, — e 2 
or -e 3 ). 

Let a G C and C a be the cycle starting from 3?a, approaching the point a parallel 
to the imaginary axis, turning anti-clockwise around a and returning to 3?a as shown 
in Figure 4. 



Im 


c 

C a 






Ra 



Re 



Figure 4- Cycle C a . 

We continue the eigenvalue E analytically along the cycle C a where a is a branching 
point which is listed in Table 2 and not listed in Table 3. The branching along the 
cycle C a is then determined as shown in Table 5. 



a = .258666 + .697448/ 
a = .224582 + .842777/ 
a = .552288 + .677536/ 
a = .314813 + .821858/ 
a = .686317+ .559106/ 


E (q) E 2 (q), E 2 (q) E (q), E 4 (q) E 4 (q), E 6 (q) => E 6 (q) 
E (q) => E 4 (q), E 2 (q) E 2 (q), E 4 (q) E (q), E 6 (q) E 6 (q) 
E (q) => E 4 (q), E 2 (q) E 2 (q), E 4 (q) £ (g), E 6 (q) => E 6 (q) 
E {q) => E 4 {q), E 2 {q) £ 2 (g), £ 4 (g) £ (g), £ 6 (g) £ 6 (g) 
E (q) E {q), E 2 {q) => E 4 {q), E 4 {q) E 2 {q), E 6 (q) => E 6 (q) 


a = .535905 + .640487/ 


Ex(q) => E 3 (q), E 3 (q) => E 1 (q), E 5 (q) E 5 (q), E 7 (q) => E 7 (q) 



Table 5. Branching along the cycle C< 



5. Convergence radius and branching points 

In section 0] we calculated the positions of the branching points of the eigenvalues 
E m {q) (m G Z) in q and described how the eigenvalues are continued along cycles. In 
this section we observe that the convergence radii of the eigenvalues E m (q) calculated 
by perturbation are compatible with the positions of the branching points. 
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For the periodic case the closest branching point from the origin is q = .258666 + 
.697448/ (|g| = .743869) and the eigenvalues E (q) and E 2 (q) are connected by con- 
tinuing analytically along the cycle C q (q = .258666 + .697448/) (see Table 5). It 
is known that the convergence radius of a complex function expanded at an origin 
is equal to the distance from the origin to the closest singular point. Hence the 
convergence radii of the eigenvalues E (q) and E 2 (q) are both .743869. 

On the other hand in section El we obtained that the convergence radii of the 
expansions of the eigenvalues E (q) and E 2 (q) around q = 0, calculated by the method 
of perturbation, are both around .749. 

Thus, convergence radii calculated by different methods are very close and com- 
patibility between the method of perturbation and the method of monodromy is con- 
firmed. Moreover, we obtain a reason why the convergence radii of the eigenvalues 
E (q) and E 2 (q) calculated in section El are very close by considering the branching 
point. To get more precise values of convergence radii calculated by perturbation, it 
is necessary to calculate more terms in k on the expansion E m (q) = E m + ^2 k Em k \ 2k 
(m = 0, 2). Generally speaking, it would be impractical to guess a convergence radius 
numerically from Taylor's expansion. 

The second closest branching point from the origin for the periodic case is q = 
.224582 + .842777/ (\q\ = .872187) and the eigenvalues E Q (q) and E A (q) are connected 
by continuing analytically along the cycle C q (q = .224582 + .842777/) (see Table 5). 
In section El we obtained that the convergence radius of the series E±(q) is around 
.875. Hence for the eigenvalue E A (q) we also obtain compatibility. 

For the anti-periodic case the closest branching point from the origin is q = 
.535905 + .640487/ (\q\ = .835115) and the eigenvalues E\(q) and E 3 (q) are con- 
nected by continuing analytically along the cycle C q (q = .535905 + .640487/) (see 
Table 5). In section El we obtained that the convergence radii of the series Ei(q) and 
E 3 (q) are both around .838. For the eigenvalues Ei(q) and E 3 (q) we see compatibility 
and we obtain a reason why the convergence radii of E\(q) and E 3 (q) calculated in 
section El are very close by considering the branching point. 

We conclude that the convergence radii of the eigenvalues E m (q) {m = 0, 1, 2, 3, 4) 
calculated by perturbation and the locations of branching points calculated by con- 
sidering the monodromy are compatible. 

We presume that all eigenvalues E m (q) (m e 2Z> ) in an + (H) (resp. all eigenvalues 
E m (q) (m G 2Z> + 1) in (Jh_(H)) are connected by analytic continuation in q. 
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Appendix A. Proof of Proposition 12.11 
To prove Prop osit ion 12 . 1 1 we review some propositions from [S], [7j. 
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Let T be the space spanned by meromorphic doubly periodic functions up to signs, 
namely 

F — -^61,63) (A-l) 

ei,63=±l 

T ei , e3 = {f(x): meromorphic \f(x + 1) = eif(x), f(x + r) = e 3 f(x)}. (A.2) 

Let V be the maximum finite-dimensional subspace in T which is invariant under the 
action of the Hamiltonian. Then it is known that dimV^ = 2n + 1 jH]. Let P(E) be 
the monic characteristic polynomial of the Hamiltonian H (see Eq. fl2.lj) ) on the space 
V, i.e. P(E) = U7=i( E ~ E i) d E i} are eigenvalues of H on V). Then the set a d {H) 
coincides with the set of zeros of P{E). From the periodicity we have o~d{H) C o~ s (H). 

Proposition A.l. Proposition 3.5] The equation 

- 4 (n(n + l)p(x) -E)-£--2 {n{n + l)p'(x))) E{x, E) = 0, 
ax 6 ax J 

has a nonzero doubly periodic solution which has the expansion 

n 

E(x,E) = Y,b J (E)p(x) n -i, (A.3) 

3=0 

where the coefficients bj(E) are polynomials in E, they do not have common divisors, 
and the polynomial b Q (E) is monic. Moreover the function S(x,E) is determined 
uniquely. 

Proposition A.2. [5, Proposition 3.7], |3 Proposition 2.6] The function 



A(x, E) = ^TW) exp J V E ^ E E ) dX (A- 4 ) 
is a solution to the differential equation M.l\) . 

It follows from Eq. (|A~4"|) that, if P(E) ^ then the functions A(x, E) and A(-x, E) 
are linearly independent (see also the proof of Lemma 3.6]) and they form a basis 
of the space of solutions to the differential equation (JT7TJ). Note that the function 
A(x, E) is also expressed as 

A^Ej^n^^w), (A.5) 

for suitably chosen A and ai {i = 1, . . . , n) (see |H §39] or [HI §23.7]). 

From the periodicity of S(x, E) and the definition of A(x, E), we have A(x + 1, E) = 
B(E)A(x, E) for some B(E). Set A s y m (x,E) = A(x, E) - (-l) n A(-x, E). Then the 
relation HA s y m (x,E) = EA sym (x, E) is obvious. 

Proposition A.3. (c.f. §4.4]; (%) If B(E) = ±1, then the function A^ m (x,E) is 
square-integrable on (0,1). 

(ii) If P(E) ^ 0, then the function A sym (x, E) is nonzero. 

(Hi) If B(E) 7^ ±1, then P(E) ^ and any nonzero solution to Eq. il.l]) is not 
square-integrable. 
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Proof, (i) Because the exponents of the differential equation at x = are — n 
and n + 1, we have the expansion A sym (x, E) = x a (co + c\X + • • • ), where cq ^ 
and (a = — n or n + 1). From the property A sym (x, i?) = (— f ) n+1 A sym (— x, E) and 
n £ Z> , we have a = n+1. Thus the function A sym (a;, E) is holomorphic at x = 0. It 
follows from the assumption 5(E) = ±1 that A sym (a; + 1, £) = ±A sym (x, E). Hence 
the function A sym (a;, E) is also holomorphic at x — 1. Since A sym (x, E) satisfies the 
differential equation ([1.1)1 . it does not have singularity on the open interval (0,1). 
Therefore A sym (x, E) is square-integrable on (0, 1). 

(ii) It follows immediately from the linear independence of the functions A(x, E) 
and A(—x, E). 

(iii) Assume that P(E) = 0. It follows from Ea. (fA~4|) that A(x,E) 2 = E(x,E). 
From the double-perioficity of the function S(x, E), we have A(x + 1, E) 2 = A(x, E) 2 . 
Hence A(x + 1, E) = ±A(x, E) and B(E) = ±1. Therefore we have P(E) ^ under 
the assumption B(E) ^ ±1. 

Assume that B(E) ^ ±1. Then -P(-E) 7^ 0, and any solution to Eq. fjl.lJ) can 
be written as a linear combination of A(x,E) and A(—x,E). The function A(x,E) 
has poles at x = and x = 1. Let f(x) = C\A(x,E) + C2A(—x,E) be a non-zero 
square-integrable eigenfunction. The function f(x) cannot have a pole at x = 
nor x = 1 for square-integrability on (0,1). If the function f(x) is holomorphic at 
x — 0, then we have Ci = — (— l) n C*2. From the periodicity we have f(x + 1) = 
dB{E)A{x,E) +C 2 B{E)- 1 A{-x,E). Hence we have C X B{E) = -(-l) n C 2 B(E)- 1 
for holomorphy of the function f(x) at x — 1. Under the assumption B(E) 7^ ±1, we 
have Ci = C2 = and it contradicts to existence of the non-zero square-integrable 
eigenfunction. □ 

From Proposition |AH1 (iii) we have ai nt (H) C a s (H). From Proposition IA. 31 (i). (ii) 
we have o~ s (H) \o~d(H) C o~i nt (H). Combining with <Td(H) C o~ s (H) we have cr s (H) = 
cr d (H) U a int (H). Therefore we obtain Proposition 12.11 (i). To prove Proposition 12.11 
(ii), it is sufficient to show the following lemma: 

Lemma A. 4. If q(= exp^i/— It)) e K and < |g| < 1, i/ien 

a d {H) n a int (iy) = 0. (A.6) 

Proof. By Proposition 13.41 it is enough to show that Cd(H) R o"n(H) = <fi. Set 

J=-^+n(n + l)p(z + r/2). (A. 7) 

Then the potential does not have poles on R. Set 



H 



(/:R-^C measurable {0 l/C^I 2 ^ < +°°' 1 (A.8) 
I /W = /(z + 2) a.e. x, J 

From the periodicity we have cr s (H) = cr-^(I) as a set. 

If g = then a basis of eigenfunctions in the space H is {exp(m7TA/— lx)} mG z and 
we have cr^I) = {7c 2 m 2 — TT 2 n(n + 1)/3| m G Z} with multiplicity. For the case q = 0, 
the set {n 2 m 2 — Tt 2 n(n J r 1)/3| m £ Z, m > ri + 1} coincides with the set o~ s (H). The 
set o~d{H) tends to the set {n 2 m 2 — 7T 2 n(n + 1)/3| m £ Z, — n < m < as q — > 0. 
We define the set cr^(if) for the case g = by o~d(H) — {^ 2 itl 2 — 7f 2 n(n + 1)/3| m £ 
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Z, — n < m < n}. Then we can check directly that a s (H) = <Jd{H) U cr int (H) for the 
case q = 0. 

By a similar discussion to Proposition 13.31 and |5, Proposition 3.3] (see also [3]), 
it follows that all eigenvalues of/(— l<g<l)on the space H can be represented 
as E m (q) (m G Z), which is real-holomorphic in q G (—1,1), E m (0) = Tr 2 m 2 — 
n 2 n(n + l)/3 and the operator / (—1 < q < 1) forms a holomorphic family of 
type (A) (for definition see [3]). From the equation a s (H) = a d {H) U a int (H) = 
(Jd(H) U o-n(H) = cth(^) an d that elements in cx^if), an(H) and 0jj(J) are all real- 
holomorphic in q (—1 < q < 1), we have <Jd{H) = {E m (q)\ m G Z, — n < m < n} 
and <Jn(H) = {E m (q)\ m G Z, m > n + 1}. Moreover we have _E m+n+ i(g) = E m (q) 
(m G Z> ) and the multiplicity of the eigenvalue E m+n+ i (q) {m G Z> ) on the space 
H is two. 

Suppose E G <Td(H) n an(H). Then is both the eigenvalue in ad(H) C Oh(J) 
(multiplicity > 1) and the eigenvalue in an(H) C o"jj(/) (multiplicity > 2) and the 
multiplicity is summed up because the operator / (— 1 < g < 1) form a holomorphic 
family of type (A). Hence the multiplicity of the eigenvalue E is no less than three. 
However that is impossible because the dimension of the solution to the second-order 
linear ordinary differential equation (/ — E)f(x) = with the boundary condition 
f(x) G H is no more than two. Thus we obtain that if — 1 < q < 1 then <Jd{H) D 
an(H) = 0. □ 



Appendix B. 

We note definitions and formulas for elliptic functions. Let uji and be complex 
numbers such that the value to^/ui is an element of the upper half plane. 

The Weierstrass p-function, the Weierstrass sigma-function and the Weierstrass 
zeta-function are defined as follows: 

p(x) = p(x\2u l7 2u 3 ) = (B.l) 



— V 

x 2 ' V (x — 2mwi — 2nu^) 2 (2mwi + 2nuj^ 2 

(m,n)£ZxZ\{(0,0)} VV d 

a(x) = a(x\2ui,2u 3 ) =x If f 1 "^ ~TT> J 



(m,n)GZxZ\{(0,0)} 
2 



X X 

■ exp h 



C(x) 



2muoi + 2nco>3 2(2771^ + 2nuj^) 2 
a\x) 



a(x) 

Setting u} 2 = — ojj, and 



ei = p(ui), r) i = (( UJ i), (i = 1,2,3). 



(B.2) 
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yields the relations 

e 1 + e 2 + e 3 = r] 1 + r]2 + r] 3 = 0, p(x) = -C'(x), (B.3) 

(p'(x)f = 4(p(x) - ei)(p(x) - e 2 )(fp(x) - e 3 ), 

C(x + 2wi) = C(x) + 2Vi, °(x + 2wi) = -a{x)e 2 ^ x+U}i \ 

p(x + 2Wi) = p(x + Ui) - 6i = , 

p(x) - a 

where = {1,2,3}. On elliptic integrals we have 

f p(t) ds 

t ~ U i = , A ( M M V (R4) 

J ei a/4(s - ei)(s - e 2 )(s - e 3 ) 

at)-m= =_, (< = 1,2,3). 

The co-p functions (i = 1, 2, 3) are defined by 

= exp(-?7ja;)(7(a; + uji)/ (a(x)a(ui)), (B.5) 

and satisfy 

Pi {xf = p(x) - a, (^ = 1,2,3). (B.6) 

Set ui = 1/2, co>3 = r/2 and g = exp^V - It). The expansions of the Weierstrass 
p function, the Weierstrass ( function and rji in the variable q are written as follows: 

p{x) = -277! + . 27 x - ^ >^ T^Sl cos 2A;7nr ' ( RT ) 
sm (7rx) f— ^ 1 — 

C(x) = 1r} X x + - — - — r + 4tt V — sin 2knx, 

tan 7rx 1 — a 
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fc=i 



kq 



2k 



2k 
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